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Abstract 

This paper studies the diffusion approximations for queueing systems with customer 
abandonment. As queueing systems arise from various applications and the abandonment 
phenomenon has been recognized to be significant, we develop the diffusion analysis for such 
systems within a general fr amework for modeling the abandonment. In light of the recent 
work bv lDai and Hel (|201(l) . our first step is to show that, under an appropriate condition, 



the diffusion-scaled cumulative number of customers who have abandoned is asymptotically 
close to a function of the diffusion-scaled queue length process. Applying this asymptotic 
relationship, the diffusion approximations for such systems are established in a wide range 
of heavy traffic regimes, including the traditional heavy traffic regime, the non-degenerate 
slowdown regime and the Halfin-Whitt regime. The key step is to develop mappings with nice 
properties that can characterize the system dynamics by connecting the primitive processes 
such as the arrival, service and abandonment processes to the processes underlying the 
systems such as the queue length or the total head count. 



1 Introduction 

Customer abandonment has been extensively studied in various queueing models, motivated by 
their applications in different service systems. Many studies focus on diffusion analysis of the 
queueing processes in heavy traffic regimes where the arrival rate and service capacity increase 
to infinity while staying roughly balanced. The asymptotic studies in such regimes indeed make 
sense since the utilization is close to one in many real service systems which are often quite large. 
This paper focuses on providing a unified approach to the diffusion analysis for the different 
scalings of the patience time distribution, and for the various heavy traffic regimes. The unified 
approach applies not only to several classical models, which have been studied before, but also 
to a few cases which have not yet been explored in the literature. The models considered in this 
paper can be generally described as queueing systems each with a single customer class and a 
single server pool, denoted by G/GI /N n + GI. Here, N n is the number of homogeneous servers 
in the server pool. We will explain in the following section how to scale N n in the general 
framework of heavy traffic regimes. 



1.1 On heavy traffic regimes 

The queueing systems with complex dynamics are usually studied in a heavy traffic regime. 
Loosely speaking, a heavy traffic regime is designed such that a significant number of events 
(customer arrivals, abandonments and service completions) can happen in a reasonable amount 
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of time. This allows the functional law of large numbers or functional central limit theorem 
(FCLT) to be used to approximate the stochastic processes underlying the queueing system. 
This kind of probabilistic characterizations can help elucidate how the systems work, and pos- 
sibly provide insights for optimal design and control of such systems. 

The general setting for a heavy traffic regime is to consider a sequence of systems indexed 
by n. Denote by A n the arrival rate to the nth system, and assume that 
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for some constant \i > 0. So the customer arrival rate increases linearly to infinity in the heavy 
traffic regime. To balance the increasing demand, the service capacity must also increase at the 
same rate. Assume that the number of servers N n = n a for some parameter a £ [0, 1] , and each 
server can process customers at rate u n = n 1_a /i. In this way, the total service rate amounts to 
n\x and just balances out the arrivals. On a more detailed level, assume that for some /3 € R, 
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This roughly means that there is an imbalance of order ^fn between the service capacity and 
the arrival rate when the magnitude of the arrival rate is of order n. We are mainly interested 
in how the stochastic process X n (t), which counts the number of customers in the system at 
time t, oscillates around N n . 

The above heavy traffic framework with parameter a was first set out by Atar ( 20121 ). in 
an effort to bridge the so-called conventional heavy traffic regime and the many-server or the 
Halfin-Whitt heavy traffic regime. The conventional heavy traffic regime is mainly for single 
serve r queues, i.e., N n = 1 or equivalently a = 0. This re gime dates back to Kingman (|1961 



19621 ). wit h the FCLT p oint of view bein g developed by llglehart and Whittl (|1970al ibl). For 



tne r ^Lii p oint or view bein g developed Dv ligienart and wmttl i|iMYUal Jq). fot 
Glvnnl (jl990h and lwhittl (|2002h . In this regime, the diffusion limit can normally be 
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characterized using a certain function of Brownian motions. The methodology and results can 
be extended t o a g eneral network s etting, but with quite substantial efforts. See, for example 
Reiman ( 1984 ) and Williams ( 19981 ). where high-dimensional Brownian motions are involved. 

Considering the other extreme by letting a = 1 yields the Halfin-Whitt regime in which 
only the number of servers increases with n whil e the service rate of eac h server is kept fixed. 
The regime was initiated in the seminal work of Halfin and Whitt ( 198ll ). where service times 
were assumed to follow the exponential distribution. Since then there has been much effort 
expen ded to solve the problem by r elaxing the assumption of exponentially distributed service 
times. Puhalskii and Reiman (2000) extended the diffusion analysis to allow the service time to 



follow phase- type distributions, and Whitt ( 20051 ) studied the same model with a finite buffer 
by assuming the service time follows a special class of service time distribution, denoted by H^-, 
which is a mixture of an exponential distribution and a unit point mass at 0. Both of these 
studies were conducted within a Markovian framework by taking advantage of the mem ory- 
less property of exp o nenti al random variables. A breakthrough was made by iReedl (|2009h and 
Puhalskii and Eeedl (j2010h . who proposed an innovative representation of the system dynam- 
ics and a regulator mapping to establish the diffusion limit in the Halfin-Whitt regime. The 
main idea is to connect the infin ite server qu eue to the many-server queue using a regulator 
mapping. In a subsequent work, iReed (j2007h described a separate approach which views the 
many- server queue from the i dle tim e perspective and leads to an equivalent diffusion limit. Re- 
cently, Kaspi and Ramanan ( 2011bl ) characterized the diffusion limit of the many-server queue 
using a stochastic partial differential equation. The approach is based on the use of a measure- 
valued process to keep track of the residual service times of the customers in the system at 
a ny time. The framework of m easure-valued processes was first developed i n an earlier work 
of Kaspi and Ramanan ( 2011al ). which analyzed the fluid limit of the system. Jelenkovic et al 
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(120041 1 demonstrated the convergence of the steady state distribution of the model with deter- 
ministic service times. Their analysis takes advantage of the fact that there is no "overtaking" 
(meaning customers depart i n the same order as they ar r ive) in many-server queues when the 
service time is deterministic. Gamarnik and Momcilovic ( 20081 ) analyzed the stationary value 
of the queue length process assuming a lattice-valued service time distribution. 

For both the conventional and Halfin-Whitt heavy traffic regimes, the arrival and total 
service rates are of order 0(n). Existing research shows that the queue length scales as (^(n 1 / 2 ), 
thus the expected delay scales as 0(n -1 / 2 ). The difference between the two regimes lies in how 
each of them achieves the total service rate of order 0(n). In the conventional one the number 
of servers is fixed at one, while the service rate of the single server scales as 0(n), resulting in 
the service time scaling as 0(n~ 1 ). On the contrary, in the Halfin-Whitt regime the number 
of servers scales as 0(n) while the service rate of an individual server is fixed at order 0(1), 
i.e., the service time scales as 0(1). The different scalings originate from different applications. 
The conventional regime often applies to manufacturing or data package transmission over the 
Internet where the service rates of machines can be quite fast. The Halfin-Whitt regime is 
typically applied to service systems where humans are involved, because it is easier to increase 
the number of servers to achieve the total service rate of order 0(n) than to make each individual 
server work faster. 

As in the work of lAtarl (I2012I ). the slowdown is defined as the ratio between the sojourn time 
and the service time experienced by a typical customer. The above discussion shows that the 
slowdown degenerates in heavy traffic t o oo and 1 in the conventional and Halfin-Whitt regimes, 
respectively. In fact, as pointed out bv lAtarl d2012l ). the degeneration occurs for any a G [0, 1] 
except a = 1/2. The heavy traffic regime with parameter a = 1/2 is therefore referred to as 
the non- degenerate slowdown (NDS) regime. The same regime, where delay and serv ice times 
remai n comparable, was also co nsidered by iMandelbaum (|2003l ). IWhittJ (|2003l ) and iGurvichl 
(120041 ). Readers are referred to lAtarl hoii ) for very interesting examples where this regime 
becomes useful. 



1.2 On customer abandonment 



On top of the basic model we can add abandonment to allow customers to leave the system 
during waiting. This is motivated by the fact that in applications outstanding orders may be 
canceled in manufacture industries, and data packets have a high probability to be dropped if 
the waiting time is too long in a transmission channel, and customers can easily hang up at a 
call center after waiting for a while. Abandonment is modeled by assuming each customer has 
a patience time, which is a random variable. A customer abandons the system on ce his waitin g 
time exceeds his patience time. The study of customer abandonment dates back to lPalml ()1937l ). 
who identified the impatient behavior of telephone switchboard customers. 

In the current literature, there are two main streams of studies on abandonment, distin- 
guished by different scalings of the patience time. The first one keeps the patience time distribu- 
tion fixed in a heavy traffic regime. This stream can be further cla ssified depending on the heavy 
traffic regime assumed. In the conventional heavy traffic regime, IWard and Glvnnl (2003) iden- 
tified the diffusion limit as a r eflected Ornstein-Uhlenbeck process for the M/M/l + M model. 
Later. IWard and Glyn n (200 5)) extended th e result to the general G/GI/1 + GI model. In the 
Halfin-Whitt reg ime, Garnett et al.l ( 2002 ) obtained the diffusion limit for the M/M/n + M 
model. Recently, Dai et alT "( 201C ) extended the diffusion analysis to a more general G / PH /n + 
GI model by applying a general continuous map to both t he fluid and diffusion-scaled pr o cesse s 
and the random-time-change theorem. At the same time, Mandelbaum and Momcilovic ( 20121 ) 
derived diffusion approx imations by connecting the G/GI /n + GI queue to the G/GI/n queue, 
building on the work of I Reed (2009), and identifying a relation between the queue length and 
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the offered waiting time processes. In the NDS regime, I Atari (|2012l ) established the diffusion 
approximation for the model with Poisson arrivals, and exponential service and patience times. 
Results in all of the above studies share the common feature that only the density of the patience 
time distribution at the origin plays a role in the diffusion limit. Intuitively, the reason is that 
the expected delay being of order 0{n~ 1 / 2 ) causes the abandonment behavior to be dominated 
by the shape of the patience time distribution between and 0{n~ 1 / 2 ). Thus, taking the limit 
n — > oo loses all the information about the patience time distribution except for the density at 
0. 

In the other stream of the literature, the patience time distribution is instead scaled by using 
the hazard rate. For the nth system, the hazard rate function is defined by 

h n (x) = h(y/nz), 



as introduced by Reed and Ward! ( 20081 ). The motivation behind this scaling is to preserve 
more information about the patience ti me distribution rather than j ust the density at a single 
point. For a more general sca ling, see iLee and Weerasinghd (|201ll ). As was pointed out by 
Zeltvn and Mandelbaum ( 20051 ) based on a statistical study of call ce nter data, the estimate 
for the hazard rate function at a single point turns out to be unstable. Reed and Ward (2008) 
obtained the diffusion approximations for both the offered waiting time process and the queue 
length process for the G/GI/1 + GI model in the conventional heavy traffic regime. Their ap- 
proach was to use a non- linear gener a lized regulator mapping to establish the weak convergence 
results. Recently, Reed and Tezcan ( 201ll ) applied the same hazard rate scaling to study the 
diffusion limit for the queue length process of the G/M/n + GI model by proving the asymptotic 
equivalence between the queue length process and the virtual waiting time process. In both stud- 
ies the stationary distribution of the limiting processes was analyzed t o obtain approximations 
for various performance measures. Extensive numerical experiments by Reed and Tezcan ( 201ll ) 
showed that the approximations involving the entire hazard rate function generally outperform 
those that only rely on the density at the origin. In particular, when the distribution density of 
the patience time changes rapidly near the origin, using only the densi ty at zero to estimate the 
drift term of the diffusion limit yields a sizable error. See Section 4 of Reed and Tezcan ( 201ll ) 
for a detailed explanation. 

Table 11.11 summarizes the existing studies on the diffusion analysis of queueing systems 
at the process level by classifying them into three heavy traffic regimes and two scalings of 
the patience time distribution. In addition to the research focusing on process level analysis, 





No Scaling 


Hazard Rate Scaling 


Conventional 


Ward and Glvnn (20031 M/M/l + M 
Ward and Glvnn (20051 G/GI/1 + GI 
Lee and Weerasinehe (20111 


Reed and Ward (20081 G/GI/1 + GI 
Lee and Weerasinehe (20111 


G/GI/1 + GI 


G/GI/1 + GI 


NDS 


Atar (20121 M/M/^Jn + M 




Halfin-Whitt 


Garnctt et al. (20021 M/Mln + M 
Dai et al. (20101 G/PH/n + GI 


Reed and Tezcan (20111 G/M/n + GI 


Mandelbaum and Momcilovic (2012) 

G/GI/n + GI 





Table 1.1: Diffusion Approximations for Systems with Abandonment 



there have been quite a few studies focusing on the steady state for such systems in heavy 
traffic. Most of them considered the case wh ere the patience time distribution is not scaled. 
For example, IZeltvn and Ma ndelbaum (200 5j) studied var ious performance measures of the 
M/M/n + GI queue based on Baccelli and Hebuterne ( 198ll ). which analyzed the M/GI/1 + GI 
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queue. iBoxma and de Waall (119941 1 developed several approximations to the probability of aban- 



donment for the M/G I/n + GI queue and perf o rmed some simulation tests. A more detailed 
survey is presented in IZeltvn and Mandelbauinl (|2005l ). which studied the steady state of the 
M/M/n + G queue. 

1.3 Our approach 

In this paper, we aim to provide a unified approach to studying customer abandonment for 
different scalings of the patience time distribution in the various heavy traffic regimes parame- 
terized by a. We first provide a general scaling of the patience time distribution F n (-) for the 
nth system in a heavy traffic regime by assuming 

y/EF n (^=) -»■ f(x), asn^oo. (1.1) 

Jn 



This asymptotic framework can unify the cases of both no scaling and hazard rate scaling. If 
there is no scaling, i.e., F n (x) = F(x) with F(0) = and F'(0) = a, then f(x) = ax. With 
hazard rate scaling, the distribution function F n (x) = 1 — exp(— J* h{^Jnt)di) for some hazard 
rate function h(-). In this case, f(x) = Jq h(t)dt. As we can see from the existing studies, 
the function f{x) given by (jl.lj) can characterize the customer abandonment in the diffusion 
limit for both cases. The intuition of (jl.ip is based on the observation that customers' waiting 
times are of order 0{n~ l l 2 ) in the class of heavy traffic regimes we are focusing on. Thus no 
matter how the patience time distribution is scaled, what is required is a meaningful limit for 
the behavior of the patience time distribution F n {-) on the interval from to <3(n -1 / 2 ). That 
is why we scale the parameter inside of F n {-) by y/n. What the function f{x)/x captures is the 
asymptotic behavior of 

F n {x/^h) - F"(0) 
x/y/n 

which turns out to be the density at the origin if there is no scaling, and a function involving the 
hazard rate if we use hazard rate scaling. In a more general sense, other scaling methods may 
be devised to emphasize different aspects of the patience time distribution depending on the 
application. It is worth pointing out that the patience time distribution is not allowed to have 
an atom at the origin, as has been required by the existing studies, meaning that all customers 
have some degree of patience. In the case where there is a non-trivial portion of the customers 
who would leave immediately upon arrival if there is a queue, the study will require some other 
techniques which are beyond the scope of this study. 

Based on the general framework (jl.ip for the scaling of the patience time distribution, we first 
identify an asymptotic sample-path relationship between the customer abandonment process 
and the queue length process in Theorem 12.11 in the ne xt section. When (jl.ip is specialized 
to the case without scaling, our result reduces to that of Dai and He ()2O10h . The relationship 



is established by connecting the abandonment process to the offered waiting time using the 
patience time distribution. The offered or the virtual waiting time can then be related to 
the queue length using a generalization of Little's law. The challenge of such study in the 
general scaling framework (jl.ip is that the queue length processes are required to be tight 



while only stochastic boundedness is needed for the case without scaling as in iDai and He 



(2010). Tightness, in particular the modulus of continuity (j2.6p . is usually difficult to verify 



thus making Theorem 12.11 less useful in the diffusion analysis. To solve this issue, we develop 
a useful tool, Corollary 12. lj which establishes the tightness of the abandonment processes only 
based on the stochastic boundedness of the queue length processes. We later design "reflection" 
mappings which take the abandonment processes together with other processes (such as the 
arrival process), which can be shown to have good properties, as the input. With these nice 
properties in hand, the queue length (or head count) process, which is the output of those 
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mappings, are thus tight. This enables us to apply Theorem 12. II in the diffusion analysis for the 
head count and queue length processes. 

The service completion process is also an element of the input to the "reflection" mappings, 
so the next step is to characterize the service completion process using the information of the 
system size. When the service time is exponentially distributed, we have a unified approach for 
all the heavy traffic regimes parameterized by a G [0, 1). When the service time distribution is 
general, we have to tailor the analysis to the specific regime where a = 1, since the generality 
of the distribution not only requires different methods, but may also lead to different results in 
different regimes. So far, we have not been able to figure out the service completion process for 
a general service time distribution when a G [0, 1). Thus, we leave it open to future research. 

The above two steps enable us to transform the balance equation for the system into an 
equation which only involves the status of the system. To characterize the diffusion limit of 
the head count process, we develop a general approach to represent the head count process as 
the mapping of an input process involving the arrival, service completion and abandonment 
processes. Our analysis reveals that the heavy traffic regimes parameterized by a £ [0,1) and 
a = 1 are intrinsically different from the mapping point of view. So we divide the analysis 
into two cases. For the case with a G [0,1), a sequence of "reflection" mappings is developed 
in Lemma l4.1i The insight from these mappings shows that diffusion-scaled total head count 
processes cannot be negative in the limit, though for each fixed nth system it can be below 
(meaning there are idling servers). Despite the limitation of analyzing the service completion 
process based on the exponential service time distribution, our "reflection" mapping approach 
still applies in general once the service completion process can be characterized for general 
service times. For the case with a = 1, we are able to obtain the diffusion limit for the system 
when the service time is generally distributed. In this case, we also develop a regul ator mappin g 
(see Lemma 15. ip , which is essentially an extension to the one in Section 4 of iReedl (|2007l ) . 
However, it is different from that one and is new to the best of our knowledge. Summarizing 
the above, our approach to the diffusion approximation is depicted in Figure 11.11 
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Figure 1.1: The approach to the diffusion approximation. 

The rest of this paper is organized as follows. We propose a unified framework for studying 
the abandonment in Section [2j A three-step procedure (Propositions I2.1H2.3P is outlined in 
Section 12.11 and the complete proofs are given in Section 12.21 In Section El we present the results 
on diffusion approximations for many-server systems with abandonment where the number of 
servers N n = n a by making some regularity assumptions. The analysis is then divided into two 
cases. Section H] analyzes the systems in the heavy traffic regimes where a G [0, 1) and Section [5] 
analyzes them in the regime where a = 1. The mapping for a = 1 is studied in the Appendices. 

To conclude this section, we introduce some notation and definitions which will be used 
throughout the paper. All random variables and processes are defined on a common probability 
space (f^J 7 , P) unless otherwise specified. We denote by Z + , R and M+ the sets of positive 



6 



integers, real numbers and nonnegative numbers, respectively. The space of RCLL (right con- 
tinuous with left limits) functions on R + taking values in R is denoted by D(R + ,R), and the 
subspace of the continuous functions in D(R+,R) is denoted by C(R+,R). For g G D(R+,R), 
g(t—) represents its left limit at t > 0, and the uniform norm of g(-) on the interval [a, b] is 
defined by 

\\g\\[a,b] = sup \g(t)\ with ||5||[o,6] abbreviated to \\g\\ b . 
te[a,b] 

For a sequence of random elements {X n ,n G Z + } taking values in a metric space, we write 
X n =4> X to denote the convergence of X n to X in di stribution. The space D(R+,R) is 
assumed to be endowed with the Skorohod J\ -topology (see lBillingslev For a probability 

distribution function G(-) defined on R + , G c {-) = 1 — G(-). For a G R, a + = max{a,0}, 
a - = max{— a,0} and [a\ is the largest integer not greater than a. We use I5 to denote the 
indicator function of set 5 C $1. 



2 The Framework for Customer Abandonment 

Consider a sequence of single station queueing systems indexed by n G Z + . Denote Q n (t) as 
the number of customers in the queue at time t, and G n (t) as the number of customers who 
have abandoned the queue by time t, in the nth system. Clearly, G n (0) = and Q n (0) is 
the number of customers who are waiting in queue at time zero. Define the diffusion-scaled 
processes Q n = {Q n {t) : t > 0} and G n = {G n (t) : t > 0} by 

\Jn yjn 

Our objective in this section is to prove an asymptotic relationship (Theorem 12. ip between Q n 
and G n under appropriate assumptions. 

Let E n {t) denote the number of arrivals by time t in the nth system, and define the diffusion- 
scaled arrival process E n = {E n {t) : t > 0} by 

g. (t) _*■(*)- **. 

Vn 

We assume that there is a sequence of constants {A n ,n G Z + } with 

lim — = u > (2.1) 

n— >oo 77, 

such that 

E n ^ E asn->oo, (2.2) 

for some process = {-E(i) : t > 0} G C(R+,R). We now introduce a sequence of i.i.d. random 
variables {7™, i G where 7™ is interpreted as the patience time of the ith arriving customer 
in the nth system. A customer waiting in the system will leave the system without receiving 
service once his patience time exhausts. The sequence {7", i G Z+} is independent of the arrival 
process E n for each n. In this paper, we assume for technical convenience that the patience 
times of the customers who are initially in the system are infinite, i.e., the initial customers in 
the queue are infinitely patient. We denote the patience time distribution as F n {-) and assume 
that for x > 0, 

v / nF n (^) -»• fix), asn^oo, (2.3) 
\/n 



where /(•) is a nondecreasing and locally Lipschitz continuous function, i.e., for any T > 0, 
there is a constant At such that for all x,y G [0, T], 

\f(x)-f(y)\<A T \x-y\. (2.4) 
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As pointed out in Section [fl there are two special cases corresponding to different approaches 
of scaling the patience time distribution: 



• No scaling. Let F n (x) = F(x) where F is a probability distribution function with F(0) = 
and F'(0) = a. In this case, f{x) = ax. 

• Hazard rate scaling. Let F n (x) = 1 — exp(— f h(y/nt)dt) for some locally bounded hazard 
rate function h(-). In this case, f(x) = Jq h(t)dt. 

In order to obtain the asymptotic relationship (Theorem 12. lj) , the key assumption is that 
the the sequence of diffusion-scaled queue length processes {Q n ,n £ Z + } is C-tight. Namely 
the sequence is stochastically bounded, i.e., for each T > 0, 



lim limsup 



sup Q n (t) > r } = 0, 

0<t<T 



and the modulus of continuity is asymptotically small, i.e., for any e > 0, 
limlimsupPj sup \Q n (s) - Q n (t)\ > e) = 0. 

<5^0 n -s>oo I. s,t£[Q,T\,\s-t\<5 > 



(2.5) 



(2.6) 



Theorem 2.1. If a sequence of G/G/N n + GI queues satisfies (|2.ip - (|2.6p . then for each T > 0, 

0, as n — ^ oo. (2-7) 



sup 

0<t<T 



G"U)-,, { f(-Q n (s))ds 
o f 1 



The statement of the result is the same as Theorem 2.1 of bai and Hel (|2O10h who studied 
the same queueing process when patience time is not scaled. In other words, the function 
f(x) = ax where a = F'(0). However, due to the general framework ()2.3p for treating patience 
time distribution, we need a stronger condition. The additional assumption we need is (|2.6|) . 
which is mainly used to deal with the nonlinearity of the function /. 



2.1 Outline of the Proof 

The proof of Theorem 12. II is based on three properties of such queueing systems, namely Propo- 
sitions [2UH2T31 which are of independent interest themselves. We will first state these properties 
and then apply them to prove Theorem l2.ll The proofs of the three propositions are postponed 
to Section [2T2l 

Following Bai and He! (|20ld ). we introduce two notions. The first one is the offered waiting 



time uf, which denotes the time that the ith arriving customer in the nth system after time has 
to wait before receiving service for each i > 1. When Q n (0) > 0, we index the initial customer in 
the queue by 0, —1, . . . , — Q n (0) + 1, with customer — Q n (0) + 1 being the first one in the queue. 
Each uf denotes the remaining waiting time of the ith customer for i = —Q n (0) + 1, • • • ,0. 
The second notion is the virtual waiting time u> n (t), which is the amount of time a hypothetical 
customer with infinite patience would have to wait before receiving service had he arrived 
at time t in the nth system. We introduce the diffusion-scaled virtual waiting time process 
uj n = {Cj n {t) : t > 0} as 

u n (t) = Vnw n (i). 

The first property of interest is the stochastic boundedness of the virtual waiting time. 



Proposition 2.1. Under assumptions ()2.ip - (|2.5p . the sequence of the scaled virtual waiting 
times {& n ,n G Z + } is also stochastically bounded. 
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The second proposition reveals an asymptotic relationship between the abandonment process 
and the offered waiting time. 

Proposition 2.2. Under assumptions (|2.1|) - (|2.5p . for each T > 0, 



sup 

0<t<T 



E n (t) 



G' 



(t) - _L y; F n 



0, as n — > oo. 



Note that neither of the above two propositions needs the modulus of continuity to vanish 
asymptotically as in (|2.6p . The next proposition establishes the relationship between the virtual 
waiting time and the queue length. For this one, condition (|2.6|) is required. 



Proposition 2.3. Under assumptions (|2.1|M|2.6p , for each T > 0, 

0, as n — > oo. 



sup 

o<t<r 



^ n {t) - Q n (t) 

Remark 2.1. By the above proposition and the triangle inequality, for any s,t £ [0,T], 
u n {t) -0J n {s) 

Thus, the C-tightness of{Q n ,n £ Z + } implies that {w n ,n G Z + } is also C-tight. 

Proof of Theorem \2.1\ According to Proposition 12. 2| it is enough to prove that as n — > oo, 



< 2 sup 

0<t<T 



0J n (t) - -Q n {t) + - Q n {t) - Q n (s) 



sup 

0<i<T 



E n (t) 



.7 = 1 



A 1 



0. 



(2i 



After adding and subtracting a new term, we have 

E n (t) t E n {t) 



.7 = 1 



/' f f(-Q n (s))ds 
o A* 



-1= £ ^ n K) - m / 



+ /i / /(u) n ( S ))(fe - /i / f(-Q n (s))ds. 

Jo A* 



f(u n (s))ds 



1 



Thus it is enough to prove that when n — > oo 

ft 



sup 

o<t<T 1 Jo 



sup 

0<t<T 



f(Co n (s))ds- f f(-Q n (s))ds 
Jo V 

E n (t) 



0, 



1 w r* 



s))ds 



(2.9) 
(2.10) 



We first prove (|2.9p . By assumption (|2.5|) and Proposition 12.11 for any e > 0, there exists 
Ti large enough such that for all large enough n, 

e 
2" 



c { sup Q n (t) > Ti) +P{ sup Gj n {t) > Ti) < 



(2.11) 



0<t<T - 1 0<t<T 

For any 5 > 0, by the local Lipschitz continuity of /(•), we have 



f(u n (s))ds 



1 f(-Q n (s))ds >S\ 

I 1 ' 



J { sup / 

<p{ sup Q n (t)>T^\+^{ sup w n (i)>Ti)+p{ sup Cj n {t) - -Q n {t) > — — ). 

^ 0<t<T > ^ 0<t<T > ^ 0<t<T A 4 1 J 
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By Proposition 12.31 the third term on the right-hand side in the above can be less than e/2 for 
all large enough n. Thus (|2.9|) is proved by (|2,lip . 

Next, we prove (pnUD . According to Lemma 3.2 of lDai and Hel (j2010h and the monotonicity 
of the distribution function F n , 



rt -, i ft -, 

I ^F n {—£j n {s-))dE n {s) < — V F n (u;™) < / ^iF n {— Cj n {s))dE 



So it is sufficient to prove the following convergence as n — > oo, 



sup 

0<t<T 



sup 

0<i<T 



1 



V^F n {—uj n {s))dE n (s)-fi f(Cb n {s))ds 



n 



V^F n (— Z n (s-))dE n (s) - n / /(tD n (s))ds 



0, 
0. 



(2.12a) 
(2.12b) 



We only pro ve ( | 2.12al) since (|2.12b|) can be proved similarly. The idea is similar to the one 
proposed bv lWard and Glvnnl (|2005l ). Define the fluid-scaled arrival process E n = {E n (t) : t > 
0} by 

n 

Condition (|2.2|) implies that as n — > oo, 

E n => e with e(i) = fxt. (2.13) 

By Remark l2.ll we have that {(oj n , E n ),n G Z + } is C-tight. So for every convergent subsequence 
indexed by n^, 

(oj Uk , £? nfc ) =>• (a), e) as n fc ->■ oo, 

for some process Co G C(R+,R). By the Skorohod representation theorem, there exists another 
probability space (Cl,J r ,¥), as well as a sequence of processes (ijo nk ,E nk ) and (u, e) defined on 
it, such that 

(w,e) = (w,e), 

and with probability one, w nfc converges to u) in D(R+,R) and E Uk converges to e in D(R+,M). 
By the continuous mapping theorem and (|2.3p . we also know that with probability one, both 
f{u nk ) and y / r^F rafc ( i y^) converge to /(£>) as n& — )• oo in D(R+,R). By Lemma 8.3 of 

1 1 v ^fc 

Dai and Pail (|1999l ). we know that with probability one, as — > oo 



sup 

0<t<T 



ri k ~F n (—^uj nk (s))dE nk (s) - (i / /(w( S ))ds 



/((& n *(a))cfe- / /(w(s))ds 



sup 

0<t<T 



0, 
0. 



As a result, with probability one, we have as — > oo, 



sup 

0<t<T 



^ F n k{ _L ^ k ( S ))dE nk (s)-fi I fU">(s)Uls 
\J n k 



(2.14) 



Since (u nk , E nk ) = (ti nk ,E nk ), we have 



^T k j F nk {-^ob nk {s))dE nk {s) -n / fu"> (*)),/.« 



n£ / F nk (u nk (s))dE nk {s)- n f(u nk (s))ds. 
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Hence (|2.14p implies that as k — > oo, 



sup 

0<t<T 



t t-t 
!T k \ F nk (uj nk (s))dE nk (s)- n f(Cj nk {s))ds 
o Jo 



0. 



Since the above convergence to zero holds for all convergent subsequences, ()2.12ap is established. 

□ 



2.2 Proofs of Propositions E3H2H3 

In order to prove the propositions, we need the following lemma. For each 5 > 0, let 

LntJ 



Lemma 2.1. // (j2.3|) holds, then for any 5 > and T > 0, 

sup |Z$(i)| =£- 0, as n —>• oo. 

0<t<T 



Proo/. Denote p n = F n and X„ 



P<^ sup \L%(t)\ > 2e 

[0<t<T 

For fc = 0, — , [nPj - 1, let 



^ h — .fri. 1 



6 i — p n . Then for any e > 0, 

v^ 1 



max L? - ) | > 2e 

l<Jfe<[nTJ \n 



[nTj 



=fe+i 



Then 



LnTj)| >e} < ±Efa(k,[nT\)] 

= — E E (^) 2 

j=fc+i 
1 LnTJ / 

^ E E ( 1 W<^}(l-2Pn)+, 
•Pn(l ~Pn)- 



=fc+l 

1 [nT\ - k 



Hence, 



mm 



X<k<\nT\ 



in w{\L?(k,\nT\)\<e\ = 1- max P ( \L?(k, \nT\)\ > s\ 

nT\-l I J KAKlnTI-1 L J 



LnTj- 

1 \nT\ - 1 . 

> 1-^^^ Pn(l-Pn)- 



(2.15) 



(2.16) 



(2.17) 



By the Ottaviani inequality (see Page 75 of Chow and Teicher ( 20031 )) and (I2.15|) . we have 

|L?(0, LnTJ) | >e} 



sup \L n s (t)\ > 2e ) < , . , 

.°^ r J min 1 < fc < LnTJ _ 1 P{|^(A ;) LnTj)| < e} 
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It directly follows from ([2"7T5]) with k = and ([2TT7]) that 



sup |L£(i)| > 2e } < 

0<t<T 



1 [nTj 



n 



Pn(l -Pn) 



i-i^(i-p n ) 



which converges to as n — > oo according to assumption (|2,3|) . 



□ 



P roo f of Proposition \2.1[ The proof is similar to the proof of P roposition 4. 4 in an earlier version 
of Dai and Hel ( 2O10l ). First, by equation (33) on page 357 of Dai and He ( 2O10l ). we have that 
for any T G (0, O n {t)), 



v 7 V / i=E»(t)+l 



This implies that 



(t + -L) - *(t) + ^ < Q" ■ (t + -L) + LP ( J- (t + -L) ) - Lp (E» (*)) 



+ v^-L n 



■ E n t + 



Then 



sup w n (t) > T } < F{ inf 

0<t<T [0<t<T 



A n r 



L n i H — — -E n {t) + 

n ) n 



E n (t) 



Lr [E n [t + 



n-F r ' 



+ LP (E n (t)) 
E n (t+ ' 



E n {t) 



< 



E n [t + — \ -E n (t) + — 
n J n 



<p{ inf 

I 0<t<T 



+ P{ sup Q n \t + —j= I > 

o<t<T V v'n J 12 J 



P| sup 



0<t<T 



r 



~ 12 / 



ls^' F "(^)-( B "( !+ ^)- fi " (4) ) a 



12 J ' 



where is given by (|2.1j) . It follows from f)2. (|2.2|) that as n — > oo 
p| inf 



r \ A n r 

+ — ) -E n (t) + — 
n I n 



0<t<T 

Lemma |2. II and f|2. 13[) implies that as n — > oo, 

r 



sup 

0<t<T 



l£lE n [t + 



LP (E n (t)) 



By ()2.3p and (|2. 13|) . we have that as n — )• oo, 



0<t<T 



sup \fn ■ F n ( —= ) • ( LWt H — — ) - L"(t) > 



12 J 



0. 



Hence, the proposition follows from assumption (|2.5p and (|2.18|> — (|2.20|l . 



(2.18) 

(2.19) 

(2.20) 
□ 
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Lemma 2.2. Under assumptions (|2.ip ~ (|2.5p . 



lim limsupP { sup ^"(wf) > T > = 0. 



Proof. According to Lemma 3.2 of Dai and He ( 20ld ) and the monotonicity of the distribution 
function F n (-), 



sup y/nF n (uj™) > T I < P J sup ^F n {-^=u n {t)) >v\. 

0<i<E n (T) I [0<t<T V n I 



(2.21) 



Thus it is enough to prove 



1 



lim lim sup P < sup ^F n (—O n (t)) > V } = 0. 

T->oo n-j-oo 0<t<T V n 



(2.22) 



Next we note 



pj sup V^F n (4=w n (t)) >rl <p(v^F n (^ri) >rj +pJ su P a> n (t)>ri 

[o<i<T V™ J I V n J (o<t<r 

For any given e > 0, by Proposition 12.11 we can choose Ti such that 



limsupP < sup Cj n {t) > Ti ) < -. 

n->oo 0<t<T 2 



Now from (|2.3|) . for the T± fixed above, 



lim lim sup P ^ ^F n (—T 1 ) > T > = 0. 



Thus for any given e > 0, there is a Tq such that when r > Tq 



1 



lim sup P { sup y/nF n {— oj n (t)) >T} <e. 

n->oo 0<*<T V n 



This completes the proof of (|2.22p . Thus the lemma is proved due to (|2.2ip . 
An immediate consequence of the above lemma is that 



E 



sup F n «) 

0<i<E n (T) 



— > 0, as n — > oo. 



(2.23) 
□ 

(2.24) 



This w ill help to prove Lemma [2.31 b elow . which is an extension of Proposition 4.2 of lDai and He 
( 20ld ). where F n (-) = F(-). The general approach of the proof is the same whether i ?n (-)'s are 
the same or vary with n. Roughly spe aking , we ne e d to u se the martingale convergence theorem 
(cf. Lemma 4.3 of Dai and He (2010) and Whitt ( 2007 )). The key condition for applying the 
theorem is (|2.24p . We thus present the result without repeating the proof. 



Lemma 2.3. Under assumptions (|2.ip - (|2.5D 

[ntj 



sup 

0<t<T 



^7= Y,^h i <<} - r 



as n — > oo, 



where h 



is a Borel measurable function such that < h(t) < 1 for all t G 
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Similar to bai et al] feoid ). we define the process 

C(t) = inf{s > : s + u n (s) > t}. 
It is clear that ( n G D (M + , R) and is nondecreasing for each n 6 Z + . 
Lemma 2.4. Under assumptions (|'2.1|) - ([2.5|) . as n —)• oo 

sup |C(t) - t| =>- 0. 

0<t<T 

Proof. By the definition of C n (^)j f° r an y * > 0, 

o<t-C n (0 <w n (C n (t)). 

Hence, 

sup \C{t) -t\< sup w n (C n (t)) < sup u n (t). 

0<t<T 0<t<T 0<t<T 

Thus the result follows from Proposition 12.11 □ 
Proof of Proposition \2.2l According to Lemma 12.31 it suffices to show that as n — > oo, 

E n {t) 



sup 

0<t<T 



* i=i 



0. 



As a customer arriving at the system before time C n (^) must have either entered service or 
abandoned the queue by time t, we have the following relationship: 

E n (("(t)-) E n (t) 

Yl 1 - - Yl 1 {7?<w?} • 

i=l i=l 

Hence it is enough to prove that as n — > oo, 

E n (t) 

sup -= V tf » <cjV . } 0. (2.25) 

o<t<r v n ^ 1 ~ 1 
u - t - 1 v i=E n (C n (t)-)+i 

Note that 

JS»(t) JS»(t) 

sup ""/= E 1 {7"<<} = sup T= E ( 1 { 7 ' l <<} _i?n ^)) 
Ki<T \M . , , , 1 1 ~ 1 ! 0<t<T s/n . mf j-i^ , , , 1 - 1 ! 



""■f I — 

n<t<T \ n 

U - I - J v i=E«(C»(t)-)+l v j = B«(^«(t)_) + i 

+ su p ~m E 

u - f - J v i=B»(C»(t)-)+l 



(2.26) 



By Lemma l2.3l the first term on the right-hand side of (|2.26[) will converge to 0. For the second 
term, 

E"(t) 

sup -= F n (u?)< sup [E n (t)-E n (C(t)-)]- sup v^K 1 ), (2-27) 

o<t<r V™ i=S n (c „ (t) _ )+1 o<*<^ o<i<B"(T) 

which weakly converges to due to (|2.13p , Lemmas 12.21 and 12.41 Thus (|2.25p holds and the 
proof is completed. □ 
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We now present some corollaries which will be used later. 
Corollary 2.1. Under assumptions (|2.ip - ([2.5p . the sequence {G n ,n £ Z + } is C-tight. 

Proof. According to Proposition ^. 21 it is enough to show the C-tightness for {-^= X^«=i F n {u)™),n £ 
Z + }. Note that for any < s < t < T, 

-= F n W)<[E n (s)-E n (t)]- sup VEF n (u?). 

V™ i=E^s)+l 0<i<EHT) 

So the C-tightness follows from the C-tightness of E n (due to (|2.13p ) and the stochastic bound- 
edness of sup 0< j<£n( T ) ^/ni 7 "^™) (due to Lemma [272]) . □ 

Define the fluid-scaled abandonment process G n = {G n (t) : t > 0} by 

n 

Corollary 2.2. Under assumptions (|2.ip - (|2,5p . as n —> oo, 

sup u n (t) => and sup G n {t) 0. 
o<t<r o<<<t 

Proof. Note that oj n {t) = Co n {t)/y/n and G n (t) = G n (t)/y/n. So the above convergence follows 
from stochastic boundedness of {u n ,n 6 Z + } (Proposition 12. ip and stochastic boundedness of 
{G n ,n G Z + } (Corollary EU). □ 

Remark 2.2. bai and He! (|2010h also prove sup w n (t) =4> for the case where f(x) = ax, see 

0<t<T 

their Proposition 4.4. 

Proof of Proposition \2.3X By the definition of oj n (t), we have 
Q n {t + u n {t)) < E n {t + u n {t)) - E n {t) 



< Q n ({t + uj n (t))-) + ( E n {t + oj n (t)) - E n {t + u n {t) - -) 

V n 

E n (t+w n (t)) 

+ Y 1 {7r<o- 

i=E n (t) 



(2.28) 



Note that 



1 \n 

-= {E n (t + Lu n (t)) - E n (t)) = E n {t + uj n (t)) - E n (t) + — • ■ uj n (t), (2.29) 

E n {t+cu"(t)) E n (t+u n (t)) 

-7= Y 1 {-/?<<} = -fa Y ( 1 {7?<<i - 

V i=B»(i) V i=B»(t) 

E«(t+w n (t)) 

+ — £ F>»). (2.30) 

v i=E n (t) 



By £22]) and Corollary E2 



sup \E n (t + u: n (t)) - E n (t)\ =>0, as n ^ oo. (2.31) 

0<t<T 
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By pTIjl and Proposition I2TTI 



— • sfn ■ co n (t) - iiu n (t)\ 0, as n ->• oo. 

n 



(2.32) 



It follows from and (|2~2j) that 



sup 

0<i<T 



£ n (i + w n (t)) - £ n (t + w n (t) - -) 



n 



< sup 

0<t<T 



E n (t + w n (t)) - E n {t + w n (i) - -) 

n 



+ 



as n — > oo. 



(2.33) 



Note that the inequality (|2.27p also holds with (C n (*) — >*) replaced by (t,t + u n (t)), so by 
Corollary 12.21 as n — ^ oo, 

E n {t+uj n {t)) 



sup 

0<t<T 



7Tr E F " W) ^°' 



(2.34) 



i=E n {t) 

Lemma [2T3l (|2.30p and (|2.34[> imply that as n — > oo, 

B"(t+w n (t)) 



sup 

0<t<T 



l 



0. 



i=E n (t) 



By condition (|2.6p . as n — > oo 



sup 

0<t<T 



Q n (t + o;™(t))-Q n ((t + a; n (t))-) 



(2.35) 



(2.36) 



Applying the above convergence f|2.31fl — (|2.36j) to the inequality (|2.28p yields that as n — > oo, 



sup 

0<i<T 



Q n (t + u n (t))-/j£j n (t) 

By condition (|2.6p and Corollary 12.21 as n — > oo, 

sup Q n {t + u n {t)) -Q n (t) 

0<t<T 

Thus, the result of this proposition follows. 



0. 



0. 



□ 



3 Diffusion Approximations 

We now consider the diffusion approximation for a sequence of G/GI /N n + GI systems indexed 
by n. For the nth system, the number of servers is 



N n = n a , 



(3.1) 



where a £ [0,1] is a parameter. Denote X n (t) as the total number of customers at time t in 
the nth system. Its diffusion-scaled process X n = {X n (t) : t > 0} is defined by 



X n (t) 



X n {t) - N n 



n 



(3.2) 



It is clear that the queue length process Q n studied in the previous section satisfies Q n = (X r 
Our main results (Theorem 13.11 (i) and (ii)) characterize the asymptotic behavior of X n . 
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In Section [2J we have established an asymptotic relationship between the queue length and 
abandonment processes with only assumptions on the arrival processes f)2. ip — (|2.2p . patience time 
distributions (|2.3|) - (|2.4|) . and the queue length processes (|2.5p - (|2.6p . Some of the intermediate 
results in Section [2] actually only require (|2.ip - (|2.5p . We will take advantage of the results in 
Section [2] to obtain the diffusion approximation for the total head count process X n . To this 
end, some assumptions on the service time are required. Considering the nth system, let vf, 
i = 1, 2, . . . be the service time of the ith arriving customer. For i = — X n (0) + 1, • • • , — Q n (0), 
vf denotes the remaining service time of the ith customer initially in service. For i = —Q n (0) + 
1, • • • ,0, vf denotes the service time of the ith customer initially waiting in queue. Customer 
—Q n (0) + 1 is the first in the queue, customer — Q n (0) + 2 is the second, and so on. We assume 
{v™ , i > —X n (0) + 1} is a sequence of independent random variables, and is independent of the 
patience times {7J 1 , i £ Z + } and the arrival process E n for each n. 

We separate the discussions of service times into two cases. Recall that is defined by ()2.1|) . 
When < a < 1, we assume that the service times follow an exponential distribution. 

Assumption 3.1. The customers' remaining service times and service times {v™, i > —X n (0)-\- 
1} are independent and exponentially distributed with rate fi n = n 1- "^. 

When a = 1, we allow the service times to follow a general distribution function H. Let H e 
denote the associated equilibrium distribution, i.e., 

H e (x) = p I (1 - H{u))du, x > 0, 
J o 

and let M denote the associated renewal function, i.e., M satisfies the following renewal equation 

M(t) = H(t) + I H(t- s)dM(s). 
Jo 

Assumption 3.2. The customer service times {vf,i > — Q n (0) + 1} are independent and 
identically distributed with distribution function H{-) which has mean /i. And the customers' 
remaining service times {t>™, —X n (0) + 1 < i < — Q n (0)} are independent and identically dis- 
tributed with distribution function H e {-). 

With Assumptions 13. ll and 13.21 the total service capacity of the nth system is n/i, regardless 
of the choice of a. We assume the following heavy traffic condition, 

— X n 

f3 n ■= 1) -». p as n -> oo, (3.3) 

n/i 

for some j3 € R. In particular, the heavy traffic condition implies (|2.1|) . We also assume the 
convergence of initial states, 

X n (0) ^( as n -»• oo, (3.4) 

for some random variable ^. 

Theorem 3.1. Assume that conditions (|2.2p - (|2.4p and (|3.3p - (|3.4p hold. For the stochastic 
processes {X n ,n £ Z + } associated with the sequence of G/GI/N n + GI systems, 

(i) if < a < 1, Assumption 13.11 holds and £ > with probability one, then X n ^ X as 
n — > oo, where X is the solution to the following 

X(t) = f + E(t) - y/jiS(t) +pnt-p f f ^ X ^ + )ds + L(t), (3.5) 

Jo " 

X(t) > 0, (3.6) 

L(t) is nondecreasing and L(0) = 0, (3-7) 

X(s)dL(s) = 0, (3.8) 
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and S is a standard Brownian motion independent of E and £; 

(ii) if a = 1 and Assumption 13.21 holds, then X n =^ X as n — > oo, where X is the solution to 
the following 

X(t) =£ + tf(t) - S(t) + + r • (A** - M (t)) 

/"* - t* (X(s)) + ( 3 - 9 ) 

+ / (X(t-s)ydM(s)-v f(^—^—)ds, 

Jo Jo V 

and S is a Gaussian process, which is independent of E and £, with the covariance given 
by 

E[S(s)S{t)] = 2 J ^M(u) -u + ^Jdu + M(s - u)M{t - v)dH{u + v) (3.10) 

for any < s < t. 

Remark 3.1. By Proposition 4-9 of \Lee and Weerasinahe hod ), we know that for any y(-) £ 



D(R+, R) and a nondecreasing and locally Lipschitz continuous function g(-) defined on R+ with 
g(0) = 0, there exist unique x(-) and £(■) in D(R + ,R) such that 

x(t) = y(t) - ! g({x(s)) + )ds + £(t), (3.11) 



o 



x{t) > 0, (3.12) 

£(■) is nondecreasing and £(0) = 0, (3.13) 

/•oo 

and / x{s)d£(s) = 0. (3.14) 



Jo 

And by Lemma [5. 1L for any a£l, y(-) £ D(R + ,R), and a nondecreasing and locally Lipschitz 
continuous function g(-) defined on R+ with g(0) = 0, there exists a unique x(-) in D(R + ,R) 
such that 

x(t) = a + y(t) + a~ ■((rt-M(t))+ [ x(t - s))~dM(s) - [ g{(x{s)) + )ds. 

Jo Jo 

Thus X in Theorem 13.11 (i) and (ii) is well-defined. 

4 G/M/N n + GI Systems with < a < 1 

Recall that G n (t) is the number of customers who have abandoned the queue before time t. 
The evolution of the process X n can be characterized by the system dynamics equation 

X n {t) = X n {0) + E n (t) - s(n n J\ xn ( s ) A N n )daj - G n {t), 

where S(-) is a Poisson process with rate one. As a result, we have 
X n {t) — N n = X n (0) -N n + E n (t) - X n t 
Si 



(/z n J (X n (s) A N n )ds) - //> J (X n (s) A N n )ds\ - G n {t) ^ 



+ (A n - nn)t + ii n \ (X n (s)- N n )-ds. 



o 
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Applying the diffusion scaling for X n , E n and G n and the definition of f3 n in (|3.3|) . we obtain 



where 
with 



X"(l) :!•"(/)+//"/ (A-"(.s))-r/.s-// / /(±(X"( S ))+)d S , 



y n (t) = X n (0) + £ n (i) - S n (t) - G n {t) + /?V, 



S n (t) 







S U" / (X n (s) AN n )ds) -iJL n / AiV n )(is 



G-(t) = G n (t)- M f f(-X n {s) + )ds. 



o 



/i 



(4.2) 
(4.3) 

(4.4) 
(4.5) 



Similarly, let 



yn = |yn (t) ; t > 0} , S n = {S n (t) : t > 0} , G n = {G n (t) :t>0}. 



To prove the convergence of process X n , for any given function g(-) defined on 
introduce a sequence of "reflection" mappings {<3?™ , n G Z + }, where for each n, <3?" : D(R 
D(R + ,R) is defined by &^(y) = x with x being a solution to 



we 



x(t)=y(t) + ti n / (x(s))-ds 



g((x(s)) + )ds. 



(4.6) 



Note that 3>™ depends on n through //• 



n 



l-a 



\x based on Assumption 13.11 The following 



lemma characterizes the sequence of such mappings and allows fi n to be more general. 

Lemma 4.1. (i) If the function g(-) is a nondecreasing and Lipschitz continuous with g(0) = 0, 
then for each n and y G D(R+,R), there exists a unique solution to (|4.6p . (ii) Let {y n ,n G Z + } 
be a sequence of stochastic processes whose paths lie almost surely in D(R+,R). Assume that 
lim n _ s>00 [i n = oo and lim infn^oo y n (0) > with probability one. If the sequence {y n ,n G Z + } 
is C -tight, then the sequence {x n ,n G Z + } given by x n = <I>™(y ra ) is also C -tight, and any limit 
(x*,y*) of a converging subsequence of{(x n ,y n ),n G Z + } satisfies (iia) x*(t) > 0; and (iib) If 
i*(t) = x*(t) + Jq g(x*(s) Jr )ds — y*(t), in other words, £* is the limit of a converging subsequence 
of {fi n f*x n (s)-ds,n£ Z+}, then J* x*(s)d£*(s) = 0. 

Proof. First we define a sequence of mappings {^,n G Z + } with ty™ : D(R + ,R) — > D(R + ,R) 
by *$>g(x)(t) = fj, n x~(t)—g(x + (t)) for each x G D(R + ,R). Clearly, \P™ is Lipschi tz continuous for 
each n since g(0) = and g(-) is Lipschitz continuous. Therefore, by Lemma 1 of lReed and Ward 
( 20041 ) . there exists a unique solution to (|4.6|) . Thus we have (i). 

Now we prove (ii). Note that by Lemma 1 of Reed and Ward (2004), if the input process 
y n G D(R + ,R), then the resulting process x n G D(R + ,R). We define z n to be the solution to 



z n (t) = y n (t)-n n / z n {s)ds 



(4.7) 

Again, Lemma 1 of R eed and Ward (l2004l l guarantees the existence and uniqueness of z n . Since 
g(-) is Lipschitz continuous, there exists no such that g{x) < pi n x for all x > and n > uq. 
Since 



x n (t) - z n (t) = fi n / (x n (s))-ds 



g((x n (s)) + )ds + /i n / z n {s)ds 
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x n {-) — z n (-) is differentiable and 
d 



dt 



x n (t) - z n (t)} = -g(x n (t)+) + fj, n [x n (ty + z n (t)] 
>-vT[x n (t)-z n (t)]. 



As a result, f t [(x n (t) - z n (t)) exp(/x n t)] > 0. Since (x n (0) - z n (0)) exp(/x"0) = 0, we have 
(x n (t) - z n (t))exp(// n i) > 0. This implies that 



x n (t) > z n (t) for all t > and n > n . 
The integral equation (|4.7|) can easily be solved to obtain 



(4.8) 



z n {t) = exp(-n n t) 
Performing integration by parts yields 



y n (0) + / eMv n s)dy n (s) 
Jo 



y n (0) + f exp(^s)dy n (s) = y n (t) + fi n f \y n (t) - y n (s)]exp(n n s)ds. 
Jo Jo 



This implies that for any < 5 < t, 



z n (t) = y n (t)exp(- f i n t) + ^exp(- f i n t) / (y n (t) - y n (s)) exp(^s)ds 

Jo 

>y n (t)exp(-fi n t)- sup \y n (t)-y n (s)\-2exv(- fJ , n 5)- sup \y n (s)\. 

t-5<s<t 0<s<t-S 



Hence, for any e > 0, 



sup (z n (t))- > e } < P{ sup (y n (t))- > 
te[o,T] I |te[o,<5] 



+ P sup (y-(t))--exp(-^t)>- 
[te[8,T\ 6 

+ P I sup sup \y n (t) -y n (s)\ > ■ 
[te[o,T] t-S<s<t 



(4.9) 



+ F{ sup 2exp(- / u n J) • sup \y n {s)\ > - } . 

[t£[0,T] 0<s<t-6 3 J 

Note that the C-tightness of {y n ,n £ Z + } implies the C-tightness of {(y n )~,n £ Z + }, and 
liminfn^oo y n (0) > implies that lim n _ 5 . oo (y n (0))~ = 0. Hence, we have 



lim lim sup P< sup (y n (t)) > - > = 0. 
<5->o n ->oo lie [0,6] 3 I 



(4.10) 



Similarly, the C-tightness of {y n ,n £ Z + } and {{y n ) , n E Z + }, and lim n ,_ >00 pO 1 = oo make the 
last three terms in (|4.9p converge to zero. Therefore, by (|4.10D . 



It follows from (Oft that 



limsupP^ sup (z n (t))- > e \ = 0. 



sup {x n (t)) =4> as n — > oo. 
te[o,T] 



(4.11) 
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Define 

%(y)(t) = x(t) and * g (y)(t) = l(t), (4.12) 
where x(-) and £(■) are given by (|3.11j) — (|3.14|) in Remark 13, 11 Now note that 

(x n (t))+ = l,"(t) + (*»)"- f g((x n (s))+)ds + ^ f V \s))~ ds . 

Jo Jo 

Hence, 

In view of (|4.1ip . we have that 

{y n + {x n )-,neZ + } isC-tight. (4.13) 
By Proposition 4.9 of Lee and Weerasinghe ( 201ll ). the mappings <3? 9 and ty g are Lipschitz 



continuous. By the continuous mapping theorem, (|4.13p implies that {(x n ) + ,n G Z + } is C- 
tight, hence {x n ,n G Z+} is C-tight by (|4.1ip . 

(iia) follows directly from (|4.1ip . Suppose (iib) does not hold, then there exists a point t 
such that x*(t) = e > and £*(■) is increasing at t. This implies that there exists a 5 > such 
that for all large enough n, f^t + 5) > £ n (t - 5) and x n (s) > e/2 for all s G (i — 5, t + S). On the 
other hand, however, when x n (s) > e/2 for all s G (t — 5, t + <5), we have that (x n ) _ (s) = for 
all s G (t — S,t + 5). This consequently implies £ n (t — 5) = £ n (i + 5). So we get a contradiction. 
Hence, (iib) holds. □ 

In view of (|4.2p and ()4.6p . it is clear that X n = $™(F n ). A key step in applying Lemma B~D 
to obtain the diffusion limit is to show the C-tightness of {Y n ,n G The involvement of the 

abandonment process in Y n makes the task challenging. To this end, we first present a study 
on the systems without abandonment. The result is not only of independent interest, but also 
paves the way for studying the systems with abandonment. 

4.1 Systems without Abandonment 

In this section, we consider the G/M/N n system. The arrival process to the nth system is still 
E n . Customers will not leave the queue until they have obtained service, so the abandonment 
process is equal to 0. To differentiate, let X^t) be the total number of customers at time t in 
the nth system without abandonment. Similarly, the process Xq = {Xq (t) : t > 0} is defined as 
Xq (t) = (Xq (t) — N n )/y/n. According to our convention given in Section [2] (see the paragraph 
between ()2.2[) and (|2.3p ). we have 

A7(0) = X n (0). 

The objective of this section is to establish the weak convergence of {Xq,ti G Z + } (see Propo- 
sition [O}. Let = {S%(t) :t>0} with 

S^(t) = -j= S(n n [\xS(s) A N n )ds) - u n f\x%{s) A N n )ds . 



10 ' JO 

In such systems, the diffusion-scaled balance equation ()4.2p becomes 

XS(t) = Y n (t) + u" [\x-(s))-ds, (4.14) 
Jo 

where 

Y n (t) = A7(0) + E n {t) - SS(t) + /3V. 

We specialize the function g = in Lemma 14.11 and denote by <3?q the corresponding reflection 
mapping for the nth system. So we have Xq = <I>Q(Y n ). 
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Proposition 4.1. Assume that conditions (j'2.'2p . (|3,3|) - (|3.4p . and Assumption \3.1\ hold. For the 

sequence ofG/M/N n systems with < a < 1, as n — > oo ; 

where So is a standard Brownian motion independent of the initial state and the arrival process. 
Proof. By (|3.ip and Assumption 13.11 



(Xq( s ) AN n )ds < nt. 



Thus 



Ft sup 

[te[o,T] 



So(*) 



> T I < P J sup -= \S{nt) - nt\ > r 
I I *e[o,T] V n 



(4.15) 



(4.16) 



which vanishes as r — > oo by the cen tral limit theorem for Poisson processes (a special case 
of Theorem 17.3 of IBillingslevI (|l999l )). So {S%,n G Z+} is stochastically bounded. Other 
components of {Y n ,n G Z + } are stochastically bounded because of the assumptions (j2.2[) . 
(I33|)- (1331). Noting that 

sup (xS(t)r < sup (y n (*))~, 

t&[0,T] te[o,T] 
we have stochastic boundedness for {sup tg [ 0j T](^o n e ^+1- Furthermore, by 



(*0 (*)) + = W) + (*?(*))" + / (*0 (*))~^ 

JO 

we have (XJ) + = <J> g (F n + (X^) - ) w ith g = (see (14 121) V F rom the Lipschitz continuity of the 
mapping & g ( ) (see Proposition 2 of iReed and Ward! $200$ )). we have stochastic boundedness 
for {(X£)+,n G Z+}. In view of = (Xfi) + - (XJ})~ , therefore, {X$,n G Z + } is also 
stochastically bounded. The stochastic boundedness immediately implies that 



where X, 



Xq =^ as n — > oo, 
i j - { XR{t) ~ Nn : i > 0}. We now prove that 

fJ, n (XQ (s))~ds =^0 as n oo. 



(4.17) 



(4.18) 



It follows from (|4.14p that 



E n {t)-\ n t (\ n -nn)t 



n 



n 



+ 



S (> J \XS(s) A iV n )d S ) - M « A JV n )d 



(4.19) 



7? 



By (|4.17j) . (I2.2p and (|3.3p . the first four terms on the right hand side of the above converge to 
zero in distribution. By (|4.15p - (|4.16p . the last term in (|4.19p also converges to in probability. 

By (X# A N n )/n = N n /n - (X^)~ , KTEh is equivalent to 

'■X$(s)AN n 



n 



-ds =>■ e(-), 



where e(-) is give n bv (12.131) , the p roposition follows from the FCLT for renewal processes (cf. 
Theorem 17. 3 of IBillingslevI ( 19991 )) and the random-time-change theorem (cf. Corollary 1 of 
Whittl <|l98Ch ). □ 
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Proposition 4.2. Assume that conditions (|2.2|) . (|3,3|) - (|3.4p . and Assumption ^. 1\ hold. For the 

sequence ofG/M/N n systems with < a < 1, 

Xq =4> Xq as n — ^ oo, 

where Xq is the solution to the following 

X (t) = £ + E(t) - So(pt) + ppit + L (t), (4.20) 

X (t) > 0, (4.21) 

Lq is nondecreasing and Lq(0) = 0, (4.22) 

X (s)dL (s) = 0, (4.23) 



o 

where Sq is given by Proposition 14.11 and is independent of E and £. 

Proof. Conditions ([2JI, ([33]) -([331), and Proposition O imply that {Y^,n G Z+} is C-tight 

and as n — > oo, 

Y n ^( + E- V/^o + 0e, (4.24) 

where e(-) is given by (j2. 13|) . By the proof of Lemma [4.11 with g = (see (|4.1ip and (|4.12p ). we 
know that 

(X»)~ =^ as n ^ oo, (4.25) 

and 

= * oW + (X T ) and f{X^{s))-ds = + «DW- 

JO 

The theorem directly follows fro m the L i pschit z continuity of $o an d ^o, the random-time- 
change theorem (cf. Corollary 1 of lWhittJ (|198d H. and (j4T24D - (l4T25l . □ 



4.2 Systems with Abandonment 

We now return to the study of systems with abandonment. Proposition 14.21 not only remains 
as an independent result, but also helps to establish a useful bound for studying s ystems with 



aband onment. Recall that Q n = {X n ) + and Q% = (Xg)+. By Theorem 2.2 of iDai and He 
(2010), on each sample path, 



Q n {t) < Q%{t) for all t > 0. (4.26) 

Since Proposition 14.21 has already shown the stochastic boundedness of {Qq , n G Z + }, by (|4.26p . 
{Q n ,n G Z + } is stochastically bounded. 

Similar to the study for systems without abandonment, the stochastic boundedness of the 
queue length processes helps to prove the following convergence for the diffusion-scaled service 
process S n defined in 



Proposition 4.3. Assume that conditions (|2.2p . (|3.3j) - (j3.4|) . and Assumption \3. l\ hold. For the 

sequence of G/M/N n + GI systems with < a < 1, as n —> oo, 

S n => y/JiS, (4.27) 
where S is a standard Brownian motion independent of the initial state and the arrival process. 
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Proof. The proof is essentially the same as that for Proposition 14.11 so we only point out the 
differences. 

Since {Q n , n G Z + } is stochastically bounded, we immediately know that (X n (-)—N n ) + /n => 
as n — > oo. Thus, in view of the fact that (X n (-) — iV n )~/n < N n /n — > as n — > oo, we have 

X n (-)-N n 

— => as n -> oo. (4.28) 



n 



Next we prove 



It follows from (fiTj) that 



■ a n (X n (s) - N n )~ , 
p v K J ^ds ^0 as n ^ oo. (4.29) 



n 



* /j, n (X n (s) — N n )~ X n (t)-N n X n (0)-N n E n (t)-X n t (X n -nfi)t 

as = 

o n n n n n 

S (V Jj X n (s) A N n ds) - fi n J * X n (s) A iV n ds 

+ — 

n 

G n (t) 



(4.30) 



+ 



The only difference between (|4,30p and (|4.19p is the term G n (t)/n, which converges to zero 
in distribution by Corollary 12.21 So (|4,29p holds by the same argument as (|4.18p . The same 
application of the FCLT and the random-time-change theorem then yields (|4.27|) . □ 

Proof of Theorem 13. II (i). The sequence of the queue length processes is stochastically bounded 
by (|4.26p and Proposition l4.21 Since /(•) is continuous, the sequence {/j, f((X n (s)) + /fj,)ds, n G 
Z + } (see ()4.5p ) is C-tight. Then the C-tightness of the sequence {Y n ,n G Z + } follows from 
condition (|3.4|) on the initial states, condition (|2.2|) on the arrival process, Proposition 14.31 for 
the service process, Corollary 12.11 for the abandonment process, and the heavy traffic condition 
(|3.3p . Note that the above reasoning relies only on stochastic boundedness of the queue length 
process, not on the modulus of continuity being asymptotically small (see Figure fl~T|) . 

Applying Lemma [4. II we have the C-tightness of {X n ,n G Z + }. This implies the C-tightness 
of {Q n ,n G By Theorem 12. 1\ we have that G n (t) =4> as n — > oo. This implies that as 

n — > oo, 

Y n £ + E - yfjiS + 0e. (4.31) 
By the proof of Lemma 14.11 with <?(•)= fif(-/n) (see (|4.11|) and (|4.12|) ). we know that 

(X n y ^0 as n ^ oo, (4.32) 

and ^ 

X n = $ 9 (r + {X n )~) and fi n [ (X n (s))-ds = ^ g (Y n + (X n )-)(t). 

Jo 

The theorem directly follows fr om the Lipsch itz continuity of & g and $> g , the random-time- 
change theorem (cf. Corollary 1 of lWhittl (|l980h ). and (l43T|) - (^21 ). □ 



5 G/GI/N n + GI Systems with a = 1 

Now we consider the case where a = 1, i.e., the number of servers = n. Denote S n (t) as 
the number of customers leaving the nth system by finishing their services. We then have the 
following simple balance equation 

X n (t) = X n (0) + E n {t) - S n (t) - G n {t). (5.1) 
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The asymptotic behavior of E n under diffusion scaling follows directly from the assumption on 
the arrival processes. In what follows, we obtain the asymptotic behavior of X n by studying 
that of S n and G n under diffusion scaling. 

The asy mptoti c anal ysis o f S n follows from the ideas of Reedl ( 20071 ). with some technical 
results from Reed ( 2009) a nd Krichagina and Puhalskiil ( 1997 ). We first present a set of rela- 
tionships from Reedl ( 20071 ) which is applicable for all G/GI/N n + G systems. Recen t ly, th e 
same set of relationships is also applied to study a multi-class model in iReed and Shakil (|2012l ) . 
Let K n (t) be the number of customers who have entered service by time t. Denote by the 
ith jump time of the counting process K n . Define 



M n (t) 



K n (i)-Q n (0) 

£ 

t=-O n (o)+i 



{ K ?+v? >t} 



and 



r(0) 



j=-X»(0)+l 

According to Propositions 2, 5 and 6 of [Reed (120071 ). the service completion process is 



(5.2) 



(5.3) 



S n (t) = S n (t) + (Xfi An) ■ fit + (X£ - n)~ ■ M(t) - / (X n (t - s) - n)~dM(s), (5.4) 



where 



S n (t) = -{W Q a (t) + M n (t)) 



(Wg(t -s) + M n (t - s)) dM{s). 



(5.5) 



This relationship was initially obtained in Sections 2 and 3 of Reedl ( 2007 ) in the context of 
a multi-server queue without abandonment by manipulating system dynamic equations in an 
ingenious manner, without involving any asymptotic analysis. It is, however, still true for queues 
with abandonment, since what (|5.4p and (|5.5p reveal is a relationship between how customers 
enter service - the process K n , and how they complete service - the process S n . The process 
K n plays the role of separating the buffer and the server pool. Of course, due to the differences 
between models with and without abandonment, the asymptotic study of K n may need different 
methods. Plugging (I5.4p into (|5.ip yields 



(5.6) 



X n (t) = X n (0) + E n (t) - G n (t) - S n (t) 

- (Xfi A n)ut - (JtJ - n)" • M(t) + [ (X n {t - s) - n) - dM(a). 

JO 

Thus (|5.6p implies that 

X n (t) = Y n (t)+ [\x n (t-s))-dM(s)-ti f f(-(X n (s))+)ds, 
Jo Jo V 



where G n is defined as in (|4.5p and 

Y n (t) 
S n (t) 



X n (0) + E n (t) - G n (t) - S n (t) + /?V + ^"(0)" • Gut - M(t)), 
S n (t) 



n 



(5.7) 

(5.8) 
(5.9) 



To further study the process X n , we introduce a reflection mapping in the following lemma. 
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Lemma 5.1. Assume that g is a locally Lipschitz continuous function with g(0) = 0. For any 
y G D(IR + ,R), there exists a unique solution x to the following equation 



x(t) = y{t) 



(z(t- s))~ dM{s)+ / g((x(s)) + )ds 



(5.10) 



Moreover, the mapping &Mg '■ D(R+,R) —> D(R + ,M) defined by x = &M,g{y) is Lipschitz con- 
tinuous in the topology of uniform convergence over bounded intervals, measurable with respect 
to the Borel a -field generated by the Skorohod J\-topology. 

The proof of this lemma will b e presented i n Appendix[Al Note that the mapping 3>M,g is an 
extension of the one in Section 4 of lReedl (|2007l ). Following this lemma, we have X n = &M,g(Y n ) 
with g(t) = fif(t/fi). The next proposition treats the convergence of diffusion-scaled service 
processes jointly with the arrival process and initial status. 



Proposition 5.1. Assume that conditions ([2.2 
sequence of G/GI/n + GI systems, 



(p£5 j) - (pr3)l and Assumption\£B hold. For the 
(X n (0),E n ,S n )^(Z,E,S), (5.11) 



where £, E and S are independent of each other, and S is a Gaussian process with continuous 
sample paths, zero mean and covariance function given by 



ES(t)S(t + 5) = 2 



{M{u)-u+-)du+ / / 
o * Jo Jo 



M(t-u)M(t + 5-v)dF(u + v) for t,5>0. 



This result is essentially the same as Proposition 9 of lReedl ((20071 ) . and in fact the proofs for 
the two are almost identical. What prevents us from directly citing the result is the definition 
of the term M n {t), which involves the process K n {t) of customers entering service. For a multi- 
server queue with or without abandonment, the process of customers entering service will be 
different. Let the processes M n = {M n {t) : t > 0} and K n = {K n {t) : t > 0} be defined 
by M n (t) = M n {t)/\/n and K n (t) = K n (t)/n, r e spect ively. According to the separation of 



variation pointed out in iKrichagina and Puhalskiil ( 19971 ). the diffusion limit of {M n ,n G Z + } 
is only affected by the fluid limit of {K n , n G It turns out that, according to the following 

lemma, the fluid limit of {K n ,n G Z + } is the same as that for queues without abandonment in 
the Halfin-Whitt regime. 

Lemma 5.2. Under conditions (|2.2p and (|3.3j) (j3.4|) . and Assumption \3.2\ we have that 

K n =>• e as n ->• oo, (5-12) 

where e given by (|2.13p . 



Proof. The limit of X n (0) being stochastically bounded implies th at (X n ( 0) — n)/n 



n — > oo. By Assumption 13.21 according to Theorems 3.2 and 3.3 of IZhand ( 20121 ). K n 
n — > oo. 



as 
e as 
□ 



To facilitate explaining why only the fluid limit of K n determines the diffu s ion li mit of S n , we 
briefly sketch the reasoning that originated from IKrichagina and Puhalskiil (|1997l ). According 
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to (|53j) . 



X"(t)-Q"(0) 

?™(0)+l 
X"(t)-Q™(0) 



- y 



i=-Q n (0)+l 



JO 



l {s+x < t} dV n (s,x) 



where 



K n (s)-Q n {Q) 

V n (s, x) = -= ( 1 {«?<*} ~ ■ 

V 71 -■_ ^n( )+l 



(5.13) 



Here —M n is defined in the same way as in iKrichagina and Puhalskiil (jl997l ) using stochastic in- 

tegral, and the fluid limit of K n (-) plays the same role as a(-) in Lemma 5.3 of lKrichagina and Puhalskii 
(1997). We t hus have the following diffusion limit of —M n , which is cited from Lemma 5.3 and 
Remark 4 of lKrichagina and Puhalskiil (|1997l ) without repeating the proofs. 



Lemma 5.3. If (|5.12p holds, then 



M n M as n ->■ oo, 



(5.14) 



where M is a Gaussian process with continuous sample paths, zero mean and covariance function 
given by 



EM(t)M(t + 5) 



H(t - u)[l - H(t + S- u)]fidu. 



(5.15) 



We can now prov e Propositio n 15.11 by connecting it to Proposition 9 of I Reed and 
some of the results of Reed (120091 ). 



Proof of Proposition 15.11 Based on Lemmas 15.21 and 15.31 we have M n =4> M as n — > oo. The 
negative sign in (|5. 14[) does not matter since M2 is Gaussian with zero mean. Now, define 



L nu % * 



H c (t ' 

v fin' 



According to Lemma 15.21 and the argument in Proposition 5.1 of I Reed (2009), we ha v e join t 
convergence (M n ,M n ) (M,M) as n — > 00. This implies Proposition 5.2 of Reed ( 20091 ) . 
which gives the joint convergence of 

(X n (0),E n ,W n ,M n ) =► (£,E,W(H e ),M) asn^oo. 



Note that here our W n and W(H e ) are the same as those of lReedl f|2007l . [2009). since these two 
terms depend only on the initial s tatus of the server pool. Thus, exactly the same approach 
from Proposition 9 of iReedl (|2007l ) applies here. Roughly speaking, the idea is to apply the 
continuous mapping to map (X n (0), E n , W n ,M n ) to (X n (0) , E n , S n ) . We omit the repeat of 
the same argument. □ 



Let Qq denote the qu eue length process of the many-server system without abandonment. It 
is proved in IReedl (12009! ) that {Q$,n £ Z+} is stochastically bounded. Again, by Theorem 2.2 
of bai and HeJ <j2O10h . with probability one, Q n (t) < Q%(t) for all t > 0. This implies that 
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{Q n , n G Z + } is stochastically bounded. We now can apply the similar argument for the system 
with q < 1 to the system with a = 1. The basic logic, which was depicted in Figure [TTT| is now 
further elaborated as follows. With the help of Corollary 12.11 the stochastic boundedness of the 
queue length process implies the C-tightness of the centered abandonment process {C n , n G Z + }. 
This, together with the standard assumptions, shows that the input process {Y n ,n G Z + } of 
the mapping <&M,g is C-tight. Utilizing the nice properties of the regulator map 3>M,gj we obtain 
the convergence of the output process {X n ,n G ^+}- In the following, we give the rigorous 
proof. 

Proof of Theorem 13.11 (ii) . Since the sequence of the queue length processes is stochastically 
bounded and /(•) is continuous, {// J* * f{{X n {s)) + / fi)ds,n G Z + } is C-tight. By Corollary 12. 1\ 
the abandonment process {G n ,n G Z + } is also C-tight. According to the definition of ()4.5p . 
{G n ,n G Z + } is C-tight. Then the C-tightness of the processes {Y n ,n G Z + } follows from 
condition (|3.4|) on the initial states, condition (|2.2|) on the arrival processes, Proposition 15. II for 
the service processes, and the heavy traffic condition (|3.3p . 

In order to prove the convergence of {X n ,n G Z + } and that the corresponding limit is given 
by (|3.9p . by the uniqueness of the solution to the mapping &M,g in Lemma I5TT1 it is enough to 
show that for every subsequence {X nk ,nk G Z + }, there exists a further subsequence converging 
to &M, g (Y) with 

Y = {Y(t) : t > 0} = + E(t) - S(t) + Put + r • (lit ~ M{t)) :t>0}. 

Since the sequence {Y n ,n G Z + } is C-tight, for the subsequence {Y nk ,rik G Z + }, there is a 
further convergent subsequence {Y n k,n' k G Z + }, along which 

Y< ^ Y °° as n' k oo, (5.16) 

in the Skorohod Ji-topology for some limit Y°° G C(R+,R). We now show that 

y°° = Y. (5.17) 



By the generalized Skorohod representation theorem (cf. Lemma C.l of Zhangj ( 20121 )). there 



exists a probability space (Cl,^, P), on which stochastic processes {Y n k,k G Z+} and Y°° arc 
defined such that 

Y< = Y< and Y°° = Y°° . (5.18) 

Moreover, with probability one 

Y< Y°° as n' k oo, 



in the space D(R + ,M) in the Skorohod Ji-topology. By Example 12.1 of iBillingslev (119991 ) 



we 



also have that Y°° G C(R+,R). Since the limit is continuous, the convergence (|5.16p also holds 
in the uniform topology. By the Lipschitz continuity of &M, g , we have that 

<S>M, g (Y n ' k ) -> QmJY 00 ) as n' fc ^oo, (5.19) 

in the space D(R + ,R) in the uniform topology, hence also in the Skorohod Ji-topology. In view 
of ()5.18p and Lemma l5.lt we have 

$MA Yn ' k ) £ $M, g (Y n ' k )=X n K 

Hence, by (f5T9l) . 

X< $M,g( Y °°) as n' fe -»-oo (5.20) 
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in the Skorohod Ji-topology. So, in view of $a/,<,(Y°°) £ C(R+,R), we know that {X n k,n' k G 
Z + } is C-tight. A direct consequence is the C-tightness of the queue length processes {Q Uk , n' k G 
Z + }. Hence, (|5,17|) follows directly from Theorem 1 2 . 1 1 and Proposition 15, 1[ Therefore, in view 
of (HTTP) and (ICTD we have 

X n 'k =^ ($> Mg (Y) as n' k oo 
in the Skorohod Ji-topology. This completes the proof. □ 



6 Conclusion 



This paper provides a unified approach to studying the basic queueing model with a single 
customer class and a single server pool in a wide range of heavy traffic regimes and with a 
generalized scaling of the patience time. The results in this paper encompass all the previous 
works summarized in Table [TTTI and extend our under stand i ng of the model to cases that have 
not yet been explored. Our results generalize those of I Atari (|20121 ) to G /M / 1 ^fri + GI, and also 
allow the hazard rate scaling on the patien ce time distribution. We also generalize the diffusion 
approximation of Reed and Tezcan ( 201ll ) to G/GI/n + GI. Due to technical challenges, the 
study of the service completion process in the regime parameterized by a G (0, 1) still remains 
open. 
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A A Regulator Mapping 



Proof of Lemma \5.1l Note that the special case where g 
Proposition 7 there). In other words, for any y € D(R + ,R 
such that 



e was proved by iReedl ([20071 ) (cf. 
there exists a unique x G D (R + , R) 



x (t)=y(t)+ (x(t-s)) dM(s). 



(A.l) 
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Define the mapping $m by x = &m(v), then <$>m is Lipschitz continuous in the topology 
of uniform convergence over bounded intervals, measurable with respect to the Borel a- field 
generated by the Skorohod J\ topology. Let A^f denote the Lipschitz constant of &m i n the 
topology of uniform convergence over the bounded interval [0,T]. 

In order to deal with the integral equation (|5.10|) . we consider the mapping * 5 given by 
z = ^ g (y) with 

z(t)=y(t)+ f g((^M(z)(s)) + )d S foryGD(R+,M). (A.2) 
J o 

Clearly, x = 'g(y)) (that is, $>M,g = $M ^g) is a solution to ()5.10p . Thus, to prove the 

lemma, it suffices to show that 

(a) the existence and uniqueness of the solution to (|A.2j) ; 

(b) fyg is Lipschitz continuous with the topology of uniform convergence over bounded inter- 
vals; 



(c) ^> g is measurable with respect to the Borel cr-field generated by the Skorohod J\ -topology. 
We focus our analysis on the bounded interval [0, T] for some T > 0. Let A^, be the Lipschitz 

4). 



constant of g on the interval [0, T]. Let 5 = 2/(3Ay Ay 



o 



Proof of (a): We first show the existence of a solution. Define no = and u n iteratively by 

u n+1 {t) = y(t) + [ g({<f> M (u n ){s)) + )ds 

for all n > 1. Then 

U n+ x(t) - Un(t) = [ | (/(!<!> .1/ (»„)(. ■*))+) - (!((<\>M(u tl _ v )(s]) + )\ds. 

Now we will show that 



o 



2 

\u n+ i - u n \\ jS < j J n J (-) n \\y\\^ s -i Ti+1 ^ s for j = 1, 2 • • • , [^Tj + 1. (A. 3) 



For j = 1, 



2 

\u n +i ~ u n \\s < A^Af,||ii„ - n n _i|| 5 x 5 < -\\u n - u n -x\\$. 



Since g(0) = 0, we have ||tti — uq\\s = \\y\\s < ||y||(|5-iTj4-i)5) as a result, 

2 2 

\\u n+1 - U n \\ S < (g) n ||y||(L5- 1 TJ+l)5 < "(gHlyllfl^Tj+l)*- 

Now assume that we have proved (|A.3p for j < k. Then for j = k + 1, 

k 

\\Un+l -Un\\(k+l)S < y^A^A^H^n - U n -i \\j & + &t ^T^\\ u n - U n -i || (k+l)S 

v- 2 „ „ 2 „ 

= > y plFn - Wn— l|b'<5 + g lFn - «n-l||(fc+l)<5 

i=i 

2-2 2 
- X] 3 J '^ n ~ 1 )' ? (3) n ~ 1 ||yll(L'5- 1 rj+l)5 + glK - «n-l||(*+l)* 
3=1 

2 2 

< A; fe+1 n fc (-) n ||y|| (L5 -i T j +1)( 5 + -||u n - u„_i||( fc+1)5 . 
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Since \\m - ito||(k+i)<5 < IMIfla-iTj+i)^ we have 



2 

|u n+ i - 1^11(^+1)5 < k k+1 C^2i k )(-) n \\y\\ {[5 -i Ti+1 ) S 

,2 
v 3 



< (* + l)^n k +\-r ||y|| (L ,-irj + i)*. 



Hence, we have proved (|A.3p . which implies 

oo oo 

^ ||"n+l - «n||r < \\u„,+ l - «n||(L<5- 1 TJ+l)5 
n=l n=l 

oo 

< ^(L^^J +l)^ lTJ+1 -^ lTJ+1 (l) n llyll(L 5 -Tj + i )5 

n=l d 

< OO. 

Thus, {u n } is a Cauchy sequence. As D(R+,R) is a Banach Space in the uniform metric, the 
sequence {u n } converges to the limit u*, which is a solution to t|A.2[) . 

The uniqueness of the solution is an immediate consequence of the Lipschitz continuity of 
ty g , which we prove next. 

Proof of (b): For any yi,y2 £ D(R+,R), the definition of S and (|A.2j) also imply that 

2 

ll*i/(lte) - * g (yi)\\8 < \\V2 - vib + -||^ 5 (2/2) - * g (yi)\\s, 
Hence, \\^ g (y2) - ^ g (yi)\\s < 3\\ll2 ~ 2/i||<5- Suppose, for i = 0, 1, . . . , k, 

WVgirn) - %(vi)\\i5 < (3i)i|y 2 - yiWis- (A.4) 

We now show that (|A,4j) holds for i = fe + 1. For any t 6 [0, (A; + 1)5], by the induction 
assumption and the definition of 5, we have 

k 

*■ 2 2 
||* 9 (2/2) - * s (yi)||t < ||j/2 - 2/1 IU + 3(3^112/2 - Vl\\t + zW^M - *g(Vi)\\t- 



i=l 



This implies that (jA.4j) holds for f = k + 1. Continuing the induction till k = [5 1 T\ yields the 
Lipschitz continuity property of *$> g . 

Proof of (c): Define 

Z(y,u){t)=y(t) + [ g({<S> M (u)( S )) + )ds. 
Jo 

First we prove the function 5 is measurable with respect to the Borel cr-field generated by the 
Skorohod Ji-topology in D 2 (R+,R) and D(R+,R). Define U(y,u)(t) = y(t) + g({u{s)) + )ds. 
It is clear that n is measurable (in fact, continuous) under Ji to pology. Since E(y,u) = 



H(y, ^m(u)) and <3?m is Ji-measurable (Proposition 7 of iReedl (|2007l )). the measurability of 
is proved. 

We know that ^n(y) = lim Z n (y,0), where Z n is iteratively defined by 

n— »oo 

Z n (y,u) = Z(y,Z n - 1 (y,u)), n = l,2,..., 
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with H°(y,-u) = u. According to Theorem 2 on page 14 of lChow and Teicherl fa00$ ). one can 
prove by induction that S n (y, 0) is a measurable function of y for each n. Since ^f g is the limit 
of 0) in the uni form metric on finite intervals, it is also the limit of H n (y, 0) in J\ topology. 

By Theorem 4.2.2 of budlev! <j2002h . we know that ^ g is a measurable function of y with respect 
to the Borel cr-field generated by the Skorohod Ji-topology. □ 
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